Simulation of Spread and Control of Lesions in Brain 
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A simulation model for the spread and control of lesions in the brain is constructed using a planar 
network (graph) representation for the Central Nervous System (CNS). The model is inspired by 
the lesion structures observed in the case of Multiple Sclerosis (MS), a chronic disease of the CNS. 
The initial lesion site is at the center of a unit square and spreads outwards based on the success 
rate in damaging edges (axons) of the network. The damaged edges send out alarm signals which, at 
appropriate intensity levels, generate programmed cell death. Depending on the extent and timing 
of the programmed cell death, the lesion may get controlled or aggravated akin to the control of wild 
fires by burning of peripheral vegetation. The parameter phase space of the model shows smooth 
transition from uncontrolled situation to controlled situation. The simulations show that the model 
is capable of generating a wide variety of lesion growth and arrest scenarios. 
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MS affects about one million people worldwide and 
causes physical and cognitive disability. There are three 
types of MS: relapsing-remitting, secondary progressive 
and primary progressive, that differ in the dynamical pat- 
terns of disease progression. There are as yet no known 
cures for MS. Patients with relapsing MS are currently 
treated with drugs that exert immunomodulatory effects 
and slow the progression of the disease; there are no ef- 
fective treatment options for the progressive forms of MS 

a. 



MS is postulated to be a cell-mediated autoimmune 
disease directed against myelin components of the CNS. 
Myelin is an electrically insulating phospholipid layer 
that surrounds the axons of many neurons. The disease 
is characterized by both inflammatory immune responses 
and neuro-degeneration. The prevailing hypothesis on 
MS pathogenesis is that auto-reactive T-lymphocytes, a 
cell type in the immune system, orchestrate a complex 
cascade of events that cause blood brain barrier disrup- 
tion and invasion of immunologically aggressive cells into 
the CNS. However, the exact causes of MS still remain 
unknown Q • The long-term goals of this research are to 
develop disease models that can be used to evaluate ther- 
apeutic strategies for this disease and, in this report, the 
specific focus is on evaluating a network model for MS le- 
sion dynamics. Literature survey indicates that network 
approaches have not been studied extensively for disease 
modeling in MS. 
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Conventional models for auto-immunity are premised 
on the occurrence of defects in the immune system that 
cause it to turn against the host tissue. A defect-free 
immune system, in this world view, purportedly only at- 
tacks pathogens, the external agents that cause illness or 
disease [3[. However, an alternative viewpoint has been 
advocated where auto-immunity is seen as the usual im- 
mune response, but directed against those components 
of the body which, in normal conditions, are inaccessible 
to the immune system d i, i, 0, H, @- For example, in 
the danger model, developed by Matzinger and collabo- 
rators [g], it is posited that stressed and injured tissues 
can mediate immune responses through the generation of 
appropriate "danger" signals. This is as opposed to the 
activation through recognition of external pathogenic cell 
types from host tissue in the conventional models. The 
concept of comprehensive immunity, developed by Nevo 
et al. 0, Q complements this alternate perspective; ex- 
perimental results supporting their idea have also been 
reported Q. The present network model is inspired by 
the alternative viewpoint. 



The key elements of the model consist of a pathological 
process that causes cellular damage and programmed cell 
death initiated through an inter-cellular signaling com- 
ponent. The programmed cell death deprives the patho- 
logical process of healthy tissue which is necessary for 
its propagation in space and time. In this, it resembles 
the action of firemen who burn peripheral vegetation to 
contain forest fires. Inter-cellular signaling is a key fea- 
ture of the model that allow pathologically damaged cells 
to propagate alarm signals and initiate programmed cell 
death. 
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II. MODEL 

An undirected, fixed radius random graph G(n, r), 
with n nodes (vertices) and radius of connectivity, r, 
is constructed to represent the CNS in this 2D network 
model. Fixed radius implies that nodes are connected 
only if they are within a distance of r. Biologically, the 
nodes of the graph can be viewed as representing cell bod- 
ies or functional units and the edges (bonds) of the graph 
can be viewed as axons or the interconnections between 
functional units. 

Let di be the degree of the ith node, i.e. the number of 
edges attached to it. The health status of each edge, at 
time t, is indicated by its "weight" , w(j, t) (j = 1, . . . , di), 
an integer number ranging from . . . w max ■ Edges with 
weight w max are fully functional or healthy units (as at 
the beginning of simulation), and those with weight zero, 
are dead. 

In the pathological process, the edges are damaged 
by lowering their weight by a single unit. However, in 
the programmed cell death process, edge weights are di- 
rectly reduced to zero. In the regeneration process, edge 
weights are raised by a unit. 

The pathological and regeneration processes are driven 
by probabilistic events wherein each edge in the affected 
region, in each time unit, has a certain probability pf (p[) 
of getting damaged (regenerated). In the general case, 
Pi (pi) i s a column vector of length w max containing 
the transition probabilities from one state of health to 
another. Probability of programmed cell death, p p , is 
independent of the health status of the edge. 

The functional or health status of the ith node is the 
sum over its edge weights, Sj(i) = Y^jLi w (ji The 
maximum possible value of Si is denoted by Si, which is 
realized when each w(j, t) — w max . 

A node damaged by the pathological process generates 
an alarm signal when the ratio of its health status to 
the fully healthy state (si(t)/Si) falls below a threshold, 
T a i. The signals received at the ith node are summed 
and propagated further when the summed signal strength 
reaches Si . The amplitude of the signal propagated along 
jth edge is equal to w(J,t). 

Programmed cell death is initiated at all the nodes 
where the propagated signals reach a threshold nf. The 
accumulated alarm signals in the region of programmed 
cell death, a circular region around the activated node 
of radius proportional to a parameter Cbf, get reset to 
zero. No additional signals are generated at these nodes, 
akin or otherwise, to the alarm signals generated in the 
pathologic process. 

The spread of the pathologic process is driven by the 
success rate in causing cellular damage. The fraction 
of edges damaged in a particular time step, among the 
total number of healthy edges visited, is computed as 
the rate of damage due to the pathologic process, Ri{t). 
The radius of the region affected by the pathologic pro- 
cess increased or decreased according to the formula, 
a x Ri(t) x ROIt=o, where ROIt=o is the radius of the 



region at the center where the initial lesion is seeded. In 
a similar fashion, the region of programmed cell death 
was computed as Cbf xROI t= o- Final results are, thus, 
invariant with respect to the initial lesion size. 



III. SIMULATION 

In the simulations reported here, a two-state model 
with w max = 1 has been employed, i.e., there are no inter- 
mediate states of health and the edges are either alive or 
dead. Additionally, the regeneration probability, f>\, was 
set to zero in order to focus exclusively on the effects of 
the interplay between the pathological and programmed 
cell death processes on lesion structure and dynamics. A 
few preliminary results using such a configuration was 
reported earlier [lOj. 

We have set n — 400 and chosen a uniform random 
distribution of points in the unit square [0, 1] x [0, 1] . The 
radius of connectivity was set to r = 0.2. All the results 
were also confirmed on a network of n = 4000, with r = 
0.06. Average degree strengths of the order of 10 are 
obtained in these configurations; degree distribution is 
gaussian. The pathological process was initiated at t = 
in a region with ROIt=o — 0.05 around the center at (0.5, 
0.5); for n = 4000, ROI t=0 = 0.015. 

A uniform probability of pathologic damage pd — 0.33 
was used, with a = 0.12. We varied r a i, Tbf and Cbf to 
identify the conditions under which the pathological pro- 
cess could be controlled by the programmed cell death. 
Larger values of Tbf indicate reduced sensitivity to the 
alarm signals, whereas a larger value of Cbf indicates 
that a larger area near the alerted node is subjected to 
programmed cell death. In the case of T a i, larger values 
indicate quicker firing of alarm signals. 



IV. RESULTS 

Fig. Q] shows the time series of damages caused to the 
system by both the pathological process and the pro- 
grammed cell death process. The first column of panels 
in the figure show the time course of instantaneous dam- 
ages to the system. The middle coumn of panels show 
the time course of the cumulative damages to the sys- 
tem. The last column of panels show the final state of 
the network at the end of the simulations. 

There are three typical scenarios which are illustrated 
in Fig. [H in the three rows from top to bottom. Figs.[T^,- 
c show a scenario where the programmed cell death is 
not of sufficient strength to significantly affect the patho- 
logical process. Note that the instantaneous damages 
from programmed cell death are hardly ever above zero. 
Also, it is seen from Fig. [TJd that the contribution of 
programmed cell death to the sum total of damages is 
insignificant. This situation occurs with a suitable com- 
bination of low T a i, high Tbf and low Cbf values. Figs.[TJi-f 
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FIG. 1: Time course of damages to the system by the patho- 
logical and programmed cell death processes. The first col- 
umn of panels show, separately, the instantaneous damages 
due to both the processes; dotted lines with asterisks indi- 
cate the damages due to the pathological process, and bold 
lines with filled circles indicate damages due to programmed 
cell death. The second column of panels show the cumulative 
damages with time. Again, the damages effected through 
both the processes have been separately shown (same sym- 
bols as earlier), as also the sum total damages to the system 
(square symbol). The last column of panels (color online) 
show the state of the network at the end of the simulations; 
the dotted lines indicate healthy edges (axons), full (blue) 
lines indicate edges damaged due to programmed cell death, 
and dark (red) lines indicate edges damaged due to the patho- 
logical process. For all the panels, Tbf — 0.5 and r a i = 0.7, 
while the Cbf values, for each row, top to bottom, are 0.2, 0.8 
and 1.5 respectively. 



shown in Fig. [2l where an averaging has been effected 
over the two sources of randomness. The phase space 
diagrams are the results of averaging over ten different 
networks, with the dynamics averaged over a thousand 
iterations. 




FIG. 2: Phase portraits of the dynamics. The y-axes in 
Figs. [2^,-b and the gray (color online) scale in Figs. (2}j-d in- 
dicate the fraction of the total number of simulations where 
the pathological process was successfully arrested. In Fig. [2^, 
T a i has been held fixed at the shown value while the different 
curves are, from left to right, for Tbf — 0.1, 0.3, 0.5, 0.7 and 
0.9. In Fig. [2b, the different curves, from right to left, are for 
T a i = 0.9, 0.7, 0.5, 0.3 and 0.1; Tbf = 0.5 is held fixed. The 
.S'-curves that result with increasing Cbf values indicate the 
smooth nature of the transition in dynamical behavior in the 
parameter space. 



shows a slightly more comples situation. In this case, pro- 
grammed cell death is clearly the dominant effect. The 
instantaneous damages caused by both the processes are 
consistently nonzero (Fig. [TJi) and the cumulative dam- 
ages (Fig. QJ) continue to grow. Nevertheless, we see 
that the damage continues to spread. In Figs. Hg-i, the 
pathological process has been well controlled. The in- 
stantaneous damages have fallen to zero in Fig. [Tg and 
the cumulative damages have levelled off even though the 
final state of the network shows that a large area of the 
network has been damaged. 

As seen from Fig.fTJ the time series is stochastic. There 
are essentially two sources of randomness in the model. 
Firstly, the pathological process is simulated by a bino- 
mial process wherein each edge visit will lead to successful 
damage if the generated random number falls below the 
value in p^ for that edge. Secondly, the random network 
itself is generated by the random distribution of the n 
points in the plane. The complete picture of the transi- 
tion from uncontrolled growth of the pathological process 
to the situation where the pathological process has been 
well arrested is seen in the parameter phase space graphs 



From Fig. [5J we see that the transition from uncon- 
trolled pathological process to arrested pathological pro- 
cess is smooth as Cbf is varied from low to high values. In 
Fig. T a i has been held fixed and the different curves, 
from left to right, are for different Tbf values, from 0.1 to 
0.9 in steps of 0.2. In Fig. [2d, Tbf has been held fixed and 
the different curves are, from right to left, for T a i values 
ranging from 0.1 to 0.9, in steps of 0.2. We shall denote 
by C£f(T a i, Tbf) the critical value of Cbf at which these 
S'-curves attain a value of 1, i.e. all instances of simula- 
tions result in the growth of pathological process being 
arrested. The combined picture in the three parameter 
space is presented in Figs. [2};-d, from two different per- 
spectives. The three different scenarios presented, from 
top to bottom, in Fig. [1] indicate, respectively, the three 
different parts, from left to right, of a typical S'-curve of 
Figs.[2k-b. 

As seen from Fig. [21 pathological process is always con- 
trolled if Cbf > Cgf. Nevertheless, the sum total damage 
to the system is not the same for all values of Cbf > C^j; 
in fact, the damage is greater, the larger the value of 
Cb f ■ Clearly, it is desirable to effect control of the patho- 
logical process with the least sum total damage to the 
system. With this in mind, average fractional damages 
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FIG. 3: Optimality in sum total damage to the system while 
effecting arrest of pathological process is shown here. With 
increasing Cbf values the minimum in sum total damage to 
the system occurs at about the same value as the critical value 
at which the fraction of simulations in which the pathological 
process is arrested attains unity (cf. Fig. [2^). This mini- 
mum occurs at higher values of C&/ with larger Tbf, and the 
minimum value also shifts upward. A similar situation oc- 
curs with different r a i values (not shown here) . Note that the 
y-axis values have been taken at t = 20. 



at different Cbf values have been plotted in Fig.[3l Three 
different curves for three different T& f values are shown in 
this figure; similar graphs can be constructed for different 
r a i values as well (not shown here; see [Tl[ ) . These aver- 
ages have been taken at t = 20 in each case. For Cbf < 
C^ values, the damages due to the pathological process 
as well as sum total damage are still growing and have 
not become stationary at t = 20; for Cbf > C^ values, 
the averages have become stationary. Nevertheless, these 
curves indicate that the least sum total damage to the 
system, with pathological process arrested, is obtained at 
Cbf = CfJ. For Cbf > C^, the programmed cell death is 
clearly effecting more damage than is necessary to arrest 
the pathological process. Since C^ depends on Tbf and 
T a i, it is not surprising to see that (cf. Fig. [2^l and Fig. [3]) 
lesser damage results when Tbf is small. 

It has already been seen (cf. Figs. [TJi-e) that arrest 
of the pathological process doe not necessarily occur if 
the damage due to the programmed cell death process is 
greater than the pathological process. What is necessary 
is that the programmed cell death process be able to 



encircle the region affected by the pathological process, 
and, furthermore, be able to create an envelope region 
of sufficient thickness to offset its likely growth factor, 
a x Ri(t) x ROIt=o- This is achieved in all instances 
of simulation when Cbf > C&/- Currently, mathematical 
analysis of this feature is being carried out to establish 
the relationship of C^J with r a ; and Tbf, and the results 
will be reported soon. 

V. CONCLUSIONS 

A physically motivated 2D network model was devel- 
oped for the CNS and employed to study the process 
of lesion formation and spread in MS. Inter-cellular sig- 
nalling of distress by the damaged cells is a key feature 
of the model which leads to programmed cell death get- 
ting activated in an attempt to arrest the lesion progress. 
The model demonstrates that the spread of the patho- 
logic process can be arrested by programmed cell death 
when the geometry of the damage inflicted by the latter 
leads to an envelope, of sufficient thickness, being cre- 
ated encircling the area of pathological process. Such an 
envelope of dead cells deprives the pathological process 
of healthy cells which can sustain its growth. The model 
shows a smooth transition, as parameters are varied, from 
the situations of run-away pathological process, through 
aggravated damage to the system caused by unsuccess- 
ful firing of programmed cell death, to the creation of 
successful envelope around the pathological process. 

The model complements the alternate viewpoint on 
auto-immunity which posits that cells and tissues sig- 
nal distress and activate the immune system. Such a 
viewpoint circumvents the need for the immune system 
to store information about likely pathogens and, also, 
makes it capable of acting in instances of cellular dam- 
age resulting from non-pathogenic causes. Further study 
of the model along with identification of the possible bi- 
ological constituents should enable comparisons with ex- 
periments and a more detailed exposition. 



Acknowledgements 

This work was carried out during a period spent at the 
Department of Pharmaceutical Sciences, State University 
of New York at Buffalo (UB), USA. Encouragement, dis- 
cussions and funding from Dr. Murali Ramanathan is 
acknowledged. Use of the compuational facilities at the 
Center for Computational Research, UB facilitated the 
work. 



[1] O. Kantarci and D. Wingerchuk, Curr. Opin. Neurol. 19 
(2006) 248; T. J. Murray, Brit. Med. J 332 (2006) 525. 

[2] M. Sospedra and R. Martin, Ann. Rev. Immunol. 23 
(2005) 683; V. H. Perry and D. C. Anthony, Pfulos. T. 
Roy. Soc. B 354 (1999) 1641. 



[3] F. M. Burnet and F. Fenner, The Production of Antibod- 
ies (MacMillan, Melbourne, 2nd ed., 1949); F. M. Burnet 
(ed.), Immunology: Readings from Scientific American 
(W. H. Freeman, San Francisco, 1976); E. Christ and 
A. I. Tauber, Biol. Philos. 15 (1999) 509. 



■5 



[4] G. J. V. Nossal, Antibodies and Immunity (Basic Books, 

New York, 2nd ed., 1978). 
[5] O. A. Smirnova and N. V. Stepanova, Biofizika 20 (1975) 

1095. 

[6] P. Matzinger, Ann. Rev. Immunol. 12 (1994) 991; 
P. Matzinger, Semin. Immunol. 10 (1998) 399. 

[7] U. Nevo, J. Kipnis, I. Golding, I. Shaked, A. U. Neu- 
mann, S. Akselrod and M. Schwarz, Trends Mol. Med. 9 
(2003) 88. 



[8] G. Moalem, R. Leibowitz-Amit, E. Yoles, F. Mor, I. R. 

Cohen and M. Schwartz, Nat. Med. 5 (1999) 49. 
[9] U. Nevo, I. Golding, A. U. Neumann, M. Schwarz and 
S. Akselrod, J. Theo. Biol. 227 (2004) 583. 
[10] T. R. Krishna Mohan, S. Sen and M. Ramanathan, To 

be published in Int. J. Mod. Phys. E, 2007. 
[11] T. R. Krishna Mohan, Unpublished results. 



